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Abstract 

We study the forms of the orbits in a symmetric configuration of a 
realistic model of the H2O molecule with particular emphasis on the pe- 
riodic orbits. We use an appropriate Poincare surface of section (PSS) 
and study the distribution of the orbits on this PSS for various energies. 
We find both ordered and chaotic orbits. The proportion of ordered or- 
bits is almost 100% for small energies, but decreases abruptly beyond a 
critical energy. When the energy exceeds the escape energy there are still 
non-escaping orbits around stable periodic orbits. We study in detail the 
forms of the various periodic orbits, and their connections, by providing 
appropriate stability and bifurcation diagrams. 
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Small molecules are quantum systems, but their classical study 
reveals their most important features. The method of the Poincare 
surface of section plays a prominent role in such a study. A difficulty 
that appears often in realistic systems is that one can not always 
choose a flat Poincare surface of section that intersects all periodic 
orbits. In this paper we use a curved Poincare surface of section in 
order to study a realistic model of the H2O molecule. We find the 
phase plots for different values of the energy of the molecule. Then 
we find all the periodic orbits of period 1 and 2 and compute their 
stability and bifurcation diagrams. 



1 Introduction 

Two dimensional Hamiltonian systems are often studied using the technique of 
the Poincare map, where one reduces the four dimensional, continuous time flow 
of the system to an associated two dimensional discrete map, by choosing an 
appropriate surface of section. Usually a flat surface of section is chosen and 
the orbits that cross this surface with a particular direction are the consequents 
of the Poincare map. 

For many systems a flat PSS is not suitable since there are orbits of these 
systems that do not cross the surface of section. In these cases another choice 
must be made. In this paper we study a symmetric model of the H2O molecule 
where we have chosen a curved PSS. This choice for the PSS permits us to study 
all the periodic orbits of the system for all the values of energy. 

The study of periodic orbits and phase plots in realistic models of simple 
molecules, such as the model we are studying in this paper, gives interest- 
ing information that can be compared with the behavior of the corresponding 
quantum system and with experiments jl], . In particular the relation between 
classical periodic orbits and quantum mechanical cigcnfunctions was empha- 
sized by several authors, starting with the work of Gutzwiller The most 
surprising result was that in many cases one simple periodic orbit determines 
the basic form of the eigenfunction [|| |], • Further details can be found by 
using more periodic orbits || and the asymptotic structures around unstable 
periodic orbits ||. 

A partial study of periodic orbits in the H 2 molecule has been done by 
Lawton and Child|l^, [fl|, Jaffe and Brumer|12 and by Kellman|l]|, who em- 



phasize the bifurcation of the local stretching modes from the normal stretching 
vibrational mode. These papers contain references to previous work on orbits 
in the H2O molecule. 

We describe the model of the H2O molecule in section 2, and use a convenient 
set of coordinates in defining the Hamiltonian. We calculate the equipotential 
surfaces for the symmetric molecule (section 3) and in section 4 we define an 
appropriate Poincare surface of section and discuss its properties. 

Then we study in detail the phase plots on such a Poincare surface of section 
and find the main islands of stability and the chaotic zones (section 5). We 
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calculate the main periodic orbits and their bifurcations and stability in section 
6. Finally we summarize our conclusions in section 7. 

2 A model for the H2O molecule 

For the study of the H2O molecule we have chosen a simple model that is well 
suited for classical calculations jL4|. The potential energy has the form 

Vhho = ■ f(R u R 2 ,R 3 ) + V$(Ri) + V^(R 2 ) + vj 1 2) H (R 3 ) + V^ ] HO {R u R 2l R 3 ) 

(1) 

where R\, R 2 and R% are the O-Hi, 0-H 2 and H4-H2 distances respectively. 
Distances are measured in Angstroms and energies in eV. The single-body term 
in equation ([!]) is given by 

Vg ] ■f(Ri,Ra,Rs) = 1-958 -i 

where a = 1.9018, p 1 =R 1 - 0.9572, p 2 = R 2 - 0.9572 and p 3 = R 3 - 1.5139. 
The two-body terms are given by the equations 

V^iR^^-D^l + am+a^ + a 3 r 2 3 ) exp {- fll r J , i = l,2 (3) 

where n = R. t - 0.9696, 01 = 4.507, a 2 = 4.884, a 3 = 3.795 and D 1 = 4.6211, 

V h ] h( R 3) = + a 4 r 3 + a 5 r 2 3 + a 6 r 3 3 ) exp {~a 4 r 3 } (4) 

where r 3 = R 3 - 0.7414, a 4 = 3.961, a 5 = 4.064, a 6 = 3.574 and D 2 = 4.7472. 
Finally the three-body term is given by 

3 

vgl (Ri,R2, R 3 ) = 0.01892 • P( Pu p 2 ,p 2 ) JJ(1 - tanh( 7lPl /2)) (5) 

i=i 

where 71 = 72 = 2.6, 73 = 1.5 and P(p\, p 2 , p 3 ) is the polynomial 

3 33 
P(pi,p2, Pa) = 1 + ClPi + Y Y G HPiPi 

i—1 i—1 j—i 

3 3 3 3 3 3 3 

CijkPiPjPk +X!EEE CtjkiPtPjPkPi (6) 

i—1 j—i k—j i—1 j—i k—j l—k 

The coefficients Cf, Cy ■, CV,-fc and are defined in the paper by Murrell and 
Carter Q. 
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Figure 1: Coordinates for the H 2 molecule. 



Instead of the original coordinates i?i , i?2 , -R3 we will use the scaled Jacobi 
coordinates (SJC) SR, Sr, 9 (figure [l]) that are defined by the equations 



1 

2a 
Sr — aR. 



SR=—^2R\ + 2Rl-Rl 



cos 6 — 



2a SR R 3 

where a is given by the relation 

mo + 2m,H 



4m 



The inverse relations are 



Sr^ 2 



(7) 



(8) 



Ri = y j + {a SR) 2 - SR Sr cos 6 
R 2 = + (a SR) 2 + SRSr cos 9 ^ 

a 

The potential function given by equations can be written in scaled 

Jacobi coordinates using equations (|]). 
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The kinetic energy of the H2O molecule has the form 

(sR2 + S^) 



T = — 
2/i 



PsR + 



Pi 



(10) 



where Psr, Psr and Pe are the conjugate momenta of SR, Sr and 9 respectively. 
The parameter /1 is defined by the relation 



/' 



m m 2 H 



mp + 2m,jj 2a 2 



(11) 



The Hamiltonian can now be written as the sum of the kinetic and potential 
energies 

H(P S r, Psr,Pe, SR, Sr, 9) = T(P SR , P Sr , Pg, SR, Sr) + V H ho{SR, Sr, 9) (12) 



3 Equipotential surfaces 

If we consider initial conditions with 

6 = tt/2 and P e = (13) 

we can easily see that 9 = and Pg = 0, meaning that the submanifold of the 
phase space determined by equations ( |l3| ) is invariant under the Hamiltonian 
flow and hence we can perform a standard type reduction in order to get a two 
degrees of freedom Hamiltonian system on this submanifold. 

The reduced Hamiltonian H is related to the original Hamiltonian through 
the relation 

H{P SR , Psr, SR, Sr) = H{P SR , P Sr , 0, SR, Sr, tt/2) = E (14) 

and its numerical value is the energy E of the molecule. This Hamiltonian 
describes the symmetric water molecule, where the distances O-Hi and O-H2 
are equal. 

The equipotential surfaces for the symmetric molecule are shown in figure 
||. The potential has two minima at SR = ±0.676065, Sr = 1.31082 where it 
takes the value Eq = —10.047. Between the two minima there is a saddle point 
at SR = 0, Sr = 1.6057 where the potential takes the value E\ = -8.993. For 
E < Ei orbits are confined inside the potential well around a minimum of the 
potential. For E > E\ the two wells join and orbits can pass from one well to 
the other. 

For values of energy greater than E 2 = —2.7892 orbits can escape to infinity 
following the two horizontal channels at the left and the right of figure ^. Escape 
through these channels corresponds to a configuration of the molecule where the 
oxygen atom has become unbounded and the two hydrogen atoms form an H2 
molecule. Although for E > E% most orbits escape there are still bounded orbits 
as one can see in figure ^(c) below. 



5 



2.5 
2.0 



1.5 
1.0 



0.5 
0.0 

-4-3-2-10 1 2 3 4 
SR 

Figure 2: Equipotential surfaces for the symmetric H2O molecule. The solid 
line is the Poincare surface of section. 



4 A Poincare surface of section 

A Poincare surface of section (PSS) should intersect almost all the orbits p5| . 
Such intersections occur naturally if the potential has a plane of symmetry. 
However this does not happen in the present case. On the other hand many 
orbits have a tendency to pass through a local minimum of the potential (al- 
though this is not always the case). Thus in order to find a suitable Poincare 
surface of section we calculated numerically for each value of SR the correspond- 
ing value of Sr where the potential function has a minimum. Then we fitted 
the numerically determined function Sr = / nU m(SR) by the analytic form 

Sr = f(SR) = af(a SR) (15) 

where 

f(x) = 0.740350 + 0.482392 exp(-3. 34209 x 2 ) + 0.631258 exp(-0. 112678 x 4 ) 

(16) 

and a is given by equation (|^). 

This curve is given in figure |^. We checked that the difference between the 
numerically determined curve Sr = fnnm(SR) and the fitted analytic form is 
very small. For our computations we used the fitted form, given by equations 
©and©. 

Contrary to the textbook approach we have not chosen SR and Psr as the 
two coordinates on the Poincare surface of section (PSS). If we had made such 
a choice we would have the following problem. Consider the orbits of figure 
0. Both orbits cross the PSS at the same point with coordinates SR, Sr with 
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Figure 3: The thick line gives the PSS. Orbits A and B cross the line at the 
same point with the same Psr but opposite Ps r - 



the same momentum Psr and they both have the same energy E. Then we 
solve equation ( O ) in order to determine Ps r , and we get two solutions ±Ps r - 
From figure || we can see that both solutions are valid, but they correspond to 
different orbits. This in turn means that on such a PSS different orbits would 
overlap and this is not permissible. 

In order to solve this problem we have chosen as coordinates on the PSS, SR 
and the component Pt of the momentum vector P = (Psr, Psr) tangential to the 



curve (15). Thus the Poincare surface of section is determined by the coordinate 
SR along the curve of figure || and the component Pt of the momentum. Some 
typical orbits on the plane (SR, Sr) are shown in figure ||. 
A variable that is canonically conjugate to Pt is 



SR 



Qt = / y/l + f'(x)dx (17) 





In fact the symplectic form in the variables Qt, Pt is 

n = (dQ t ) a (dp t ) (is) 

and it can be shown that this is preserved. The form (jl^) in the variables SR, 
P t is written 

n = (l + f'(SR))^ 2 {dSR) A(cLPt) (19) 

Thus the two-form dSRAdP t is not preserved in general. In our study we use the 
variable SR instead of Q t because it is not practical to compute Qt numerically. 

Although the coordinates SR and P t are not canonically conjugate, this does 
not cause any problems in the study of the system. In fact the integration is 
performed in the canonical coordinates SR, Sr, Psr, Psr and the coordinates 
SR and Pt are used only on the PSS. The result is that the phase plots on the 
PSS are only slightly distorted relative to the phase plots that we would get if 
we had chosen canonically conjugate coordinates on the PSS. 
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Figure 4: The orbits a and c belong in the islands of la and lc respectively. 
The orbit b lies between these two islands. The thick line represents the PSS. 



5 Phase plots on the Poincare surface of section 

In figures we give the distribution of the orbits on the Poincare surface of 
section for various values of the energy E. 

The energy of figure |^(a) is close to the minimum of the potential. In this 
plot practically all the orbits are ordered, defining closed invariant curves. Many 
invariant curves form two sets around two stable periodic orbits. The upper one 
is the orbit la and the lower one the orbit lc (see next section). However there 
are also invariant curves between these two sets, that start on the left or upper 
side of the boundary, and terminate on the right side of the boundary. In figure 
|] we can see why this is happening. In this figure we see three orbits for this 
value of the energy in configuration space. The orbits labeled as 'a' and 'c' 
belong to the islands of the periodic orbits la and lc respectively. We see that 
the PSS crosses the outline of these orbits at two non-adjacent sides. When 
this is happening we get closed invariant curves on the PSS. But for the orbit 
labeled as 'b', that lies between the two islands of figure |(a), the PSS crosses 
the outline of this orbit at two adjacent sides (figure |J), and on the PSS we get 
an open invariant curve. 

Figures ||(b) , 0(c) and |B|(d) refer to larger energies and are given in a larger 
scale than figure |q(a). The orbit lc is continued as the double period orbit 2c in 
figure ||(b), the second point being near the left limit of the figure. The islands 
around both points 2c are now larger. In this case also the chaotic orbits are 
insignificant. 

In figure ||(c) the orbit 2c is unstable (and is represented by two points at 
the centers of two dark regions) and it has produced by bifurcation four stable 
points that correspond to a stable orbit 4c of period 4, surrounded by islands 
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;ure 5: Phase plots for small values of the energy, (a) E = —9.9 (b) E 
E = -9.1 (d) E = -8.99. 



of stability above and below each point of the orbit 2c. In this case there are 
some small chaotic zones around each unstable orbit. 

The chaotic zones have grown considerably in figure |(d). In this case there 
are still 4 islands of stability 4c, as in figure |(c), but very small. All the 
rest of the area previously occupied by bifurcations of the orbit 2c and their 
corresponding islands are now chaotic. This case is just above the limiting 
value of the energy E\ = —8.993 at which the two potential wells join into one. 
Thus there is one symmetric figure for SR < that joins the figure ||(d) at 
the small throat near the point (0,0). The point (0,0) represents the orbit lb, 
which is unstable and produces a large chaotic region around it. 

In figure ||(a) we give the phase plot for a larger value of the energy (E = —8) 
when the two potential wells are well connected at SR = (see the equipotentials 
of figure ||). Figure ^(a) is symmetric with respect to the axis SR = 0. The 
apparent asymmetries in the islands of this figure are due to the use of non- 
symmetric initial conditions on both sides of the axis SR = 0. In this case the 
orbit la is unstable and has generated by bifurcation the orbit 2a, surrounded 
by two islands of stability on the right half of figure ^(a). A symmetric orbit la' 
and its bifurcations 2a' appear on the left side of the figure. The orbits la, la' 
have generated two chaotic domains, that surround the islands 2a (respectively 
2a'). At the lower left and upper right sides of this figure there are two elongated 
islands that belong to the family 2b. The center of the figure (0, 0) represents 
again the unstable orbit lb, that generates a large chaotic domain. 

Most of figure ^(a) is covered by one chaotic orbit for which we have calcu- 
lated 6-10 4 points. We see that the points of the chaotic orbit are not distributed 
uniformly on the PSS. This is a transient phenomenon and is due to stickiness 
of the orbits near cantori surrounding the islands 2a, or the islands 2a'. These 
cantori form partial barriers. Orbits need a long time in order to cross the 
cantori and for considerable times the regions inside and outside the cantori 
seem separated. This causes the apparent non uniformity of the distribution of 
points on the PSS. But after a larger number of iterations the distribution of 
the points becomes uniform. 

The phase plots for larger values of E have some similarity with figure 0(a) 
but there are also differences. For example in figure ^(b) we see that there are 
two symmetric islands, one on each side of the axis SR = 0. These islands 
belong to the irregular family 2b, which is symmetric with respect to the point 
(0,0) (see next section), and not to two different families like the families la 
and la' of figure ||(a). The chaotic domain in this case is again very large. 

For E above the energy E = E2 orbits can escape. In figure ^(c) we can 
see the islands around an orbit that belongs in the family 2b. The orbits that 
are outside these islands escape and this is why the area between the islands is 
empty. 

The percentage of the area of the PSS covered by organized orbits as a 
function of the energy is presented in figure 0. We see that for values of energy 
close to Ei = —10.047 (the minimum value of the potential) almost all orbits 
are organized. An abrupt change happens at values of energy around E = —9. 
For E > —9 most of the PSS is filled by chaotic orbits. There are islands of 
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Figure 7: Proportion of regular orbits, as a function of the energy E 



stability for even larger values of the energy, and we notice that for values of 
the energy close to E = — 3 the percentage of the organized orbits increases, 
although it remains below 0.2 (or 20%). This increase is due to the increase of 
the size of the island around the family 2b. The islands continue to increase in 
size even when E becomes larger than the escape energy E% = —2.79, but for 
still larger E they become smaller. 

6 Periodic orbits 

We have computed the most important families of periodic orbits for the sym- 
metric water molecule. The stability diagram of the period- 1 families as a func- 
tion of the energy of the molecule can be seen in figure ||(a). Such a diagram 
gives the Henon stability parameter [jl6| a of each family as a function of E. For 
a period-fc orbit of an area-preserving map cr, the Henon stability parameter is 
defined by the relation 

a = ^ce(D[a k }(x Q )) (20) 

where xq is a point of the periodic orbit and D[a k ](xo) is the Jacobian matrix of 
a o ■ ■ • o a (k times) evaluated at xq. When the stability parameter of a periodic 
orbit lies in the interval from —1 to 1 the orbit is stable, otherwise it is unstable. 

Since the Poincare map is not symplcctic, as we have explained in section 4, 
the determinant of the Jacobian matrix Da(x) is not equal to 1 for arbitrary 
points x on the PSS. However for a periodic orbit with period k the determinant 
of -D[cr fe ] is exactly equal to 1. In fact, using equation ( |l9| ) it can be proven that 
the determinant of the Jacobian matrix of a k at a point xq = (SRq, Pto) is given 
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(c) (d) 

Figure 8: Stability diagrams: (a) Families la, lb, lc (b) Families la, lc, 2c 
(continuation of lc) (c) Family la and its period 2 bifurcations (d) Families 2b, 
2r, 2z. 
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Figure 9: Periodic orbits la and lc for E = —10.04664 in configuration space. 
The line (a) represents the numerically determined minima of the potential and 
the line (b) represents the selected PSS. 



by the relation 



where Xk — cr k (xo) — (SRk,Ptk) is the fc-th image of xq. If an orbit is periodic 
with period k, then SR = SRk hence det D[a k ](xo) = 1. 

For values of energy close to the minimum value of the potential (E = E\ = 
— 10.047) there exist two period-1 families. We call these families la and lc. 
They correspond to the two normal modes of the water molecule. In figure |i] 
we can see the corresponding periodic orbits at E — —10.04664 in configuration 
space. 

The family lc has a peculiar behaviour that is the result of our choice of the 
PSS. This family is born at E = -10.047, but beyond the energy E = -9.331 
it appears as a period-2 family, that we call 2c (figure ||(b)). We note that no 
bifurcation appears here. In order to understand what is happening, we should 
see figure [l(]. In this figure we see one orbit of the family lc that intersects the 
PSS only at one point, as expected for small energies. An orbit of the same 
family for a larger value of the energy intersects the PSS at two points and thus 
it appears as a period-2 orbit. We will call the change of the apparent period 
of a family that is not related to a bifurcation, a transition. We note here that 
there would be a transition even if we had chosen, as the PSS, the numerically 
determined curve of the minima of the potential, instead of the fitted function. 
Family 2c becomes unstable at E = —9.164 and generates by bifurcation two 
stable period 4 families. 




(21) 
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Figure 10: The transition lc-^2c. The orbit 'a' (E = —9.8) belongs to the 
family lc. The orbit 'b' (E = —9.2) belongs to the family 2c. 



For values of the energy between E = E\ — —10.047 and E = —7.2 the 
largest islands on the PSS are around the family la and the family 2a that 
bifurcates from la (figures ^]and ^|(a)). The family la is created at the minimum 
energy of the system E = Ei = —10.047 as a stable family and becomes unstable 
at E = —8.95 (figure ||(a)). At this value of E it generates by bifurcation the 
period-2 family 2a (figure ||(c)). Notice that when the stability parameter a of 
a period- 1 orbit is equal to —1 the same orbit described twice has a stability 
parameter +1. The family 2a starts at E = —8.95 with stability parameter 
a = 1, and exists for larger values of E. The family la remains unstable until 
E = —2.47, where the period-2 family 2e is generated by bifurcation from the 
family la, existing for smaller values of E (figure ||(c)). 

The family 2a for E = —6.82 becomes unstable by crossing the line a = — 1 
and becomes stable again for E = —5.62. The stability parameter of the family 
2a takes the value a = 1 at E = —5.568. At this energy an unstable period-2 
family, that we call 2f is created. The family 2e becomes unstable (going from 
larger to smaller values of the energy) for E = —2.987 and becomes stable again 
for E = —5.661. The stability parameter of the family 2e takes the value a = 1 
at E = —5.675. At this energy the unstable family 2f, that was created at a 
larger energy from 2a, joins 2e (figure ||(c)). The bifurcation diagram of the 
families la, 2a, 2e, 2f and lb is shown in figure m This A gure gives the value 
of SR for each energy. We see that the family 2f exists only in a small interval 
of values of E and joins the family la at its maximum E and the family 2e at 
its minimum E. This explains why the families 2a, 2e and 2f form a loop in 
figure |(c). 

Between E = —2.47 and E = —2.00 the family la remains stable. At 
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Figure 11: Bifurcation diagram of family la 



E = —2.00 the family la joins the family lb, and for larger values of energy the 
family la does not exist. We may thus say that the family la is a bifurcation 
of the family lb. The family lb is represented by the point (0, 0) on the surface 
of section. It is created at the energy E^ — —8.993 (figure |](a)), at the saddle 
point of the potential (figure ||), with coordinates SR = 0, Sr = 1.6057. For 
energies somewhat larger than E2 the family lb is unstable (a — 9.73). As E 
increases, the stability parameter of the family lb becomes larger, but later (for 
larger E) it becomes smaller. At E — —2.00 this family becomes stable. It 
becomes again unstable at E = —1.50 with stability parameter a = — 1. For 
larger E the family lb is unstable, but at E — —0.0205 it becomes stable again 
until E = —0.0110, where it becomes unstable, by crossing the line a = 1 and 
thus generating another period- 1 family. As the energy approaches the energy 
E = the stability parameter of the family lb tends to infinity, and the return 
time (period) of the periodic orbit also tends to infinity. 

In figure [l2| we can see that periodic orbits of the family lb appear as straight 
lines in configuration space with SR = and varying Sr. In physical space this 
corresponds to a colincar configuration of the molecule with the oxygen atom 
in the middle and the two hydrogen atoms oscillating symmetrically at each 
side of the oxygen atom. For large intervals of the energy this configuration is 
unstable. In figure |l2| we can also see how the periodic orbits of the family la 
appear in configuration space for different values of the energy. We note that 
there is also a family la' that is symmetric to the family la with respect to 
the axis SR = 0. As the energy increases the families la and la' approach the 
family lb, and join lb for E = —2.00. 

At energy E — —8.06 two irregular period-2 families are created (figure 
||(d)). Irregular families are created in pairs (one stable and one unstable) by 
a tangent bifurcation. This means that irregular orbits are not created by the 
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Figure 12: Periodic orbits of the families lb (axis SR — 0) and la. The dotted 
line represents the PSS. 



bifurcation of an existing periodic orbit, but from each other. The stable family 
is called 2b and the unstable one is called 2r. The stability parameter of the 
family 2r grows from a = 1 to a = 17.12 as the energy increases and then 
decreases, until it reaches the value a — 2.77 at E = —6.454. For this energy 
the family 2r appears as a period-4 family and we have a transition, as in the 
case of the transition lc— >2c. 

The stable family 2b plays an important role in the structure of the phase 
space. The stability diagram of the family 2b and the period-2 families that 
bifurcate from 2b can be seen in figure ||(d). The diagram is very complex, 
much more than the stability diagram of the family la. Here we consider only 
the family 2b and its equal period bifurcations. The stability parameter of the 
family 2b crosses (or just reaches) the horizontal line a — 1 eight times. Each 
time this happens a family of period-2 orbits bifurcates from the family 2b. For 
E = 0.27 the family 2b joins a new unstable family 2z and they both disappear 
there, i.e. they do not exist for larger E. 

The families 2b, 2r and 2z are symmetric with respect to the point (0, 0). The 
families that bifurcate from 2b are not symmetric but they come in symmetric 
pairs. In figure |](d) some families that have bifurcated from the family 2b 
suddenly terminate because they are going through transitions. Some orbits of 
the family 2b can be seen in figure [l^. In physical space these orbits correspond 
to oscillations of the oxygen atom between the two symmetric minima of the 
potential. 

The island around the orbits of the family 2b is the largest island for values 
of energy close to E = —3. Generally, from E = —8 to E = —3 the island 
around 2b occupies a significant measure of the PSS. 

In the present paper we studied the stability (and instability) of the various 
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Figure 13: Orbits of the family 2b in configuration space. The dotted line 
represents the PSS. 



4 




.2 I 1 1 1 1 

-9.0 -8.5 -8.0 -7.5 -7.0 

E 

Figure 14: Stability diagrams of the family 2a for deviations on the symmetry 
plane (a) and perpendicular to it (ag). 
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types of orbits on the plane of symmetry 9 = ir/2 and Pq = 0. However orbits 
that are stable on this plane may be unstable in a direction perpendicular to 
this plane. An example is shown in figure [l4|, where we have calculated two 
stability parameters for the family 2a. One is the usual stability parameter 
a, as in figure ^(c), and the other is the stability parameter ag for deviations 
out of the symmetry plane. We see that while the family 2a is stable with 
respect to deviation on the symmetry plane it becomes unstable with respect to 
perpendicular deviations for E = —8. At the transition to instability there is a 
bifurcation of a new period 2 family, which extends out of the symmetry plane 
and is initially stable. 

The study of the stability of orbits in a direction out of the symmetry plane 
is important in many cases. Such studies in galactic dynamics have been made 
since a long timeJlT]]. A detailed study of the stability and the bifurcations of 
periodic orbits with respect to the third dimension in the H2O molecule will be 
given in a future paper. 

7 Summary and conclusions 

We have made a systematic study of the classical orbits for the symmetric 
configuration of a realistic model of the H2O molecule. 

In order to study the distribution of the ordered and the chaotic orbits for 
this system, we had to choose an appropriate Poincare surface of section (PSS). 
The usual choice for the PSS is a plane surface in phase space, but for the 
particular system under study here, this would not be convenient, because for 
any choice of a plane PSS there would be important periodic orbits that would 
not cross it. 

The PSS we chose is a surface passing approximately through the minima of 
the potential. Most non-escaping orbits cross this PSS. Choosing appropriate 
coordinates on this PSS, we were able to reduce the study of the dynamics of 
the symmetric H2O molecule to the study of a two-dimensional map. With 
this particular choice of the PSS we ensure that the Poincare map describes the 
complete dynamics of the molecule. A minor problem with this particular choice 
for the PSS is that some periodic families appear to change period without a 
bifurcation, due to an extra intersection of a periodic orbit with the PSS. 

We found the distribution of the orbits on this PSS, distinguishing between 
ordered and chaotic orbits. The ordered orbits are represented by isolated points 
if they are periodic and by invariant curves if they are quasi-periodic. Most or- 
bits are ordered for small energies. But as the energy increases beyond a critical 
value the proportion of ordered orbits decreases abruptly. When the energy in- 
creases beyond the escape energy (E2 = —2.7892), most orbits escape to infinity, 
but there are still orbits trapped around stable periodic orbits. This remains 
true even when the energy increases above E3 = 0, but for positive energies the 
regions of ordered motion are very small, and practically insignificant. 

We studied the main periodic orbits of period 1 and 2 for negative values of 
the energy. We constructed stability diagrams and bifurcation diagrams for the 
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various families of periodic orbits. These diagrams are essential in understanding 
the role of various stable orbits in trapping non-periodic orbits around them. In 
fact for different energies the trapping takes place around different orbits and we 
can separate the interval between the critical energy and the escape energy into 
two regions. In the first region (lower energies) the ordered orbits are trapped 
mainly around the families 2a and 2a'. In the second region, close to the escape 
energy the ordered orbits are trapped mainly around the irregular family 2b. 
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